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We solve the wave equation for the electromagnetic field tensors associated with the 
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O I CMBR photons in a universe with scalar metric perturbations. We show that the 

^ ■ coupling of the electromagnetic fields with the curvature associated with the scalar 

^ ■ perturbations gives rise to an optical rotation of the microwave background photons. 
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CN . The magnitude of the gravitationally generated V-Stokes parameter anisotropy Ay, is 
O 

Q I however very small compared to the linear polarisation caused by Thomson scattering. 
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In the standard treatment of the microwave background anisotropics [1] it is assumed 
that the gravitational potential at the last scattering surface gives rise to temperature 
anisotropics at large angular scales through the Sachs- Wolfe effect [2]. Thomson 
scatterings near the surface of last scattering cause a hnear polarisation of the CMB 
photons. 

The anisotropics of CMBR are governed by the the following set of equations: The 
Boltzmann equation for the photon ditribution function f{pi,PQ, z) 

DX ^ dx^ ^ dx^ dX dp' ^ ' 

the geodesic equation of photons in the perturbed metric 

^ + T%yp^ = (2) 
and the photon disperison relation 

gi.uP^p" = 0. (3) 

The dispersion relation (3) is derived in the zeroeth order of eikonal approximation 
from the electromagnetic field wave equation in a curved spacetime. In the standard 
eikonal approximation of photon trajectories in gravitational fields, the red shift of 
the photon is independent of the polarization of its spin and gravity therefore makes 
no contribution to the polarization of the CMB radiation fields. 



It has been shown in [3] that the gauge invariant wave equation of electromagnetic 
fields in curved spacetimes is given by 

V^V^F.A + Rp^.x + R\F,p - R\Fxp = (4) 

Electromagnetic field tensors couple to the Riemann and Ricci curvature tensors of 
the background spacetime and these couplings can give rise to polarisation of elec- 
tromagnetic waves in gravitational fields. In this we show that scalar perturbations 
to the Robertson- Walker metric near the surface of last scattering can give rise to 



circular polarization to the CMB radiation. We calculate the anisotropy in this grav- 
itationally generated circular polarisation as a function of angular separation of the 
beams. We find however that the magnitude of this circular polarization anisotropy 
is unlikely to be observed in the planned CMB polarisation experiments [4]. 



We consider the scalar perturbations in the Friedmann- Robertson- Walker Universe 
described by the metric [5] 

ds^ = 0^(77) [(1 + 2^p)dri^ - (1 - 2(t))'-iijdx'dx^] (5) 

In the synchronous gauge the two scalar perturbation arc the Newtonian potential ip 
and the spatial curvature (j). The eikonal equation obeyed by the photon wave vector 
% — {Qo,Q) which follows from the wave equation (1) is given by [3] 

• + = (6) 

where foi is the amplitude of the electric field associated with the radiation and the 
anisotropic gravitational permeability tensor e'^- is given by [3] 

= {R^o + R'o S) + i-^R'"o^ + ^R'' o^ f + R' r " R' o '£) (7) 

The components of the permeability tensor e'j which follows from the metric (2) are 
given in Appendix A. 

Consider a photon wavenumber q and choose the z-axis along q. The wave equations 

— * 

for the transverse E field components is 



where 

e\ = -^p(l-2V^) + ^(fcifci)^(fc,r7) (9) 

el = ^k'^k^il^ik.r]) (10) 

and similarly for the others, where k — (/ci, A;2, ^3) is the wavenumber of the gradients 
of the scalar perturbation ip^k,-!]). The dispersion relations of the two propagating 



modes E± are obtained by setting the determinant of the matrix operator (5) to zero 
to give 



The eigenvectors E± are given in terms of the mixing angle x 
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where the mixing angle x is given by 

tan 2x '- 

Expressing e*^ in terms of polar coordinates of , we have 



(13) 
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k sin 9 ip.i/2 



qo± = q[l + ^p{l~2^)T ^2.2 



(14) 



and the mixing angle as defined in (p!3[) is 



tan 2x = tan 2(j) 



(15) 



(16) 



The normal mode solutions E± are therefore given by 

E+ = [Eicos0 + E2sin0]e^((''o+^)-''^) 

E_ = [-Eisin0 + E2cos0]e*((««-'')-'?^) 
The Stokes parameters of a radiation field are described by the 2x2 complex matrix 
[5], 

(5^_1/ I + Q U-^V\ 
I 2[ U + iV I-Q ) 

where 



(18) 



where the time average is taken over duration larger than the inverse frequency. In 
unpolarised radiation Q = U = V = The degree of polarisation is defined by 
a vector P with magnitude P = a/Q^TT7^"h^T^. The quantity i = i/Q^~+T/^ 



I = [E+E^* + E^E * 
Q = (e+E* - E_E * 



U= (E+E_* + E_E^*/ , 
V = %{ E.E * - E*E_ 



denotes the degree of linear polarization and V denotes the circular polarisation of 
the radiation. Using the eigenmodes in the gravitational field given by (|T^) and using 
the definitions ( [T8| ) we find the Stokes parameters of the propagating beam are 
Qiv) = Q iVi) cos 2(j) + U {r]i)sin2(j) 

Uirj) = U {rji) cos 20 cos5{rj — rji) — Q (rji) sin 2(j) cos(Aco'(?7 — r/j)) 

(19) 

^iv) = {-U iVi) COS 2(j) + Q {r]i) sin 20) sin(Acj(77 - r]i)) 
nv) = I iVi) 

where 

Auj = (go+ - go-) = ^ sin^ t]) (20) 
is the rate at which the plane of polarisation is rotated. 

We see that if at some initial time r]i there is non-zero Q {rji) or U [rji) then at a 
later time a circular polarization V{ri) is generated. The intensity of the beam I {7]) 
remains unchanged. One can express (16) as differential equations for evolution of 
Q, f/, V as 



f- = (Auj) sin (Acur]) [U cos 20 - Q sin 20] 

i = (21) 



= (Auj) cos {Aut]) [-U cos 2(f) + Q sin 20] 



drt 
dV 
drj 

To the leading order in (Auj) we see that U = Q = and 



(9y 

— = (Ao^) [[/cos20-gsin20] (22) 

Assuming that at the time of decoupling rji there is a non-zero U (rji) Q {rji) due to 
Thomson scattering, one can estimate the degree of circular polarisiation V in CMB 
radiation. 



From the form of equation (0) it is clear that the circular polarisation V averaged 
over all angles of the gravitational perturbation vector k will be zero. The rms value 
of V defined as Av{q, k, cosO) =< (V— < V >Y >^/^. Similarly defining Aq and, 
Au) , we see that the polarization anisotropics obey the coupled set of Boltzmann 
equations [6], 

Au + ikfiAu = -kAU (23) 



/^Q-\{l-P2{li))Sp 



(24) 



Av + ifc/iAv = ( Aq sin 20) ( Acj) (25) 

where /i = cos 6' and Sp = — AT2 the quadrupole temperature anisotropy. The 
gravitational contribution arises as a source term for the V^-polarisation mode while 
the Q and U modes are generated by Thomson scattering parametrised by k = 
(Xe'iT'eO'TO'iri) / a (rjo)) with Xe the ionised fraction, rif. the electron number density and 
(Jt the Thomson scattering cross section. In the tight coupling regime (keeping only 
leading order terms in ) the solutions of ( [23| ) and ( ^4]) are given by 

Af; = (26) 
AQ = -^sin2MT, (27) 

where we see that the Q polarization is generated by a quadrupole temperature 
anisotropy A^j. 

Using ( P5| ) and ( P?]) the solution of Ay {k, r], /i) is 

Ay (fc, r/, /i) = e-^'^^''-''') f ^) sin^ ^At^ Acu (r/ - r],) (28) 



One can estimate that the two point correlations 

C''^{e)^{AT{hi)Av{n2)) (29) 

and 

C^^{e)^{AT{ni)AT{n2)) (30) 

have relative magnitudes 

^ A^ . (ry - r^,) ^ ^//-i (31) 

Taking the photons wavelength ~ Im and the smallest measurable metric perturba- 
tion to be the size of galaxies ~ lOOkpc, we find that the GTR contribution to po- 
larisation anisotrpy is smaller than the temperature anisotropy by a factor of 10~^^. 
This may be compared with the corresponding factor for the ratio of the polarization 
anisotropy due to Thomson scattering to temperature anisotropy which is 10"''. 



This ratio is very small unless the anisotropies are observed at very small angles (large 
Krnax) but here the observations are difficult because of contamination from point 
sources. In conclusion, we note that gravitational couplings of the scalar perturbations 
on photon polarisations can generate circular polarisation anisotropies in principle. 
However, in practice, it would be difficult to observe. 



Appendix 



For the metric (2) with 



tjj Riemann and Ricci components of curvature are 



ke 



i^\S\ - rk5\ + i^\5\ - i^\5\ 



I (1 - 2^) - 2^ 
\-5 



Xi X3 _ Xi S.] 



ok 



R\ 



R° 



R\ 



^{^l-3[(!-S) (1-2^)-^-^] 



and the field equations R'^ = SttG {T) - \5'jT) give 



\a a-^ 



:i-2v^)-2M_^_ 

(f + f^)(l-2^)-^-^ + ^^ 
wherein the equation of state p = ^ is used. 



AT^GpbO? (1 + 3i/) (1 + 
477(7 (pba^) (1 + I/) 

471^^602 (1 - z/) (1 + 



The e matrix is then given by 
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(^ = 1,2,3) 
no summation 
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